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Abstract 



A semi-microscopic model for the low-energy photodisintegration of the 9 Be 
nucleus is constructed, and the experimental data are analyzed with its help. 
The older radioactive isotope data are supported by this analysis. The the- 
oretical photodisintegration cross section is derived. The astrophysical rates 
for the reaction a + a + n — > 9 Be+7 and the reverse photodisintegration of 
9 Be are calculated. The new reaction rate for a + a + n — > 9 Be+7 is compared 
with previous estimations. 



PACS numbers: 25.20.Dc, 25.40.Lw, 97.10Cv, 21.60.Gx 21.45. +v, 
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I. INTRODUCTION 



Recently fully microscopic calculations of nuclei with A < 9 have become feasible JTJ^]. 
The 9 Be nucleus is such a system of special interest, as it allows tests of theories of interaction 
of composite particles |§. The properties of low-energy continuum of 9 Be are of particular 
importance in this connection. However, the corresponding experimental data on the low- 
energy photodisintegration of 9 Be are not in mutual agreement (see Fig. 1). In the present 
work we develop a semi-microscopic model to describe the process, and we analyze the 
experimental data with its help. The model accounts simultaneously for both resonant and 
non-resonant contributions to the cross section. An estimation of the reliability of various 
data is obtained and a theoretical photodisintegration cross section is derived. 

We also calculate the reaction rates of the reaction 9 Be+7 — > a + a + n and the reverse 
reaction for astrophysical conditions. These reaction rates are of relevance in the high- 
entropy bubble in type II supernovae, an astrophysical site that has been suggested for the 
r-process The baryonic matter in this bubble is dominated in the beginning by a- 

particles, neutrons, and protons. The abundance distribution shifts then to higher masses 
through the recombination of the free a-particles, neutrons, and protons. This generates the 
so-called a-process leading to the formation of massive isotopes (A « 100). The reaction 
path in the a-process is mainly determined by requirements of nuclear statistical equilibrium 
and depends also on the reaction rates of the various recombination paths bridging the mass 
5 and 8 gaps. It has been shown that there are three principal reaction paths from 4 He to 
12 C: 

(i) 4 He(2a, 7 ) 12 C 

(ii) 4 He(a n, 7 ) 9 Be(a,n) 12 C 

(hi) 4 He(2n, 7 ) 6 He(a,n) 9 Be(a,n) 12 C. 

It was shown in Refs. that the triple-alpha process (i) can be neglected compared 
to the reaction sequence (ii) via 9 Be under r-process conditions in the a-process. Also the 
reaction path (iii) via 6 He can be neglected for this scenario MM- This is true even if the 
reaction rate of 4 He(2n,7) 6 He is strongly enhanced ||, because then 6 He is also destroyed 
very effectively through photodissociation. Therefore, for the a- and r-process the reaction 
4 He(a n,7) 9 Be plays a key role in bridging the unstable mass gaps at A = 5 and A = 8. 

The reaction rates of 4 He(a n,7) 9 Be and the reverse photodisintegration of 9 Be were 
estimated in Ref. || from the experimental photodisintegration cross section. However, 
Ref. || did not include information on which experimental data their estimate was based. 
In view of the astrophysical relevance of these reactions we recalculate in the present work 
the rates of the first step of the reaction (ii) above. 

The same problem is also addressed in Ref. || . These authors obtain the resonant contri- 
bution to the 9 Be(7,n) 8 Be cross section from the Breit- Wiener formula for the first excited 
state of 9 Be with the parameters taken from Ref. |TD[ .[] In order to calculate the non-resonant 



contribution they introduce a single-particle potential with the depth chosen to reproduce 
the ground state, calculate both ground- and final-state continuum wave functions in this 



1 We note that the T and Eq parameters of the resonance used in Ref. |J seem to be incorrect. 
The resonant properties of the l/2 + state of 9 Be will be considered in our future work. 
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potential, and multiply the cross section obtained by the shell-model spectroscopic factor. 
They then add this cross section constructively or destructively to the resonant cross sec- 
tion to establish possible upper and lower bounds for the reaction rates. This procedure 
has certain shortcomings: a resonant contribution to the cross section should not emerge as 
an addition to the dynamic model used, since a correct quantum mechanical model should 
necessarily contain such a contribution itself, along with the non-resonant contribution and 
an interference term. Besides, the potential wells used for the ground state and continuum 
state should in fact be different: an additional spin-orbit potential, for example, should be 
present in the ground p-state as compared to the continuum s-state. 

In our model we use a three-body specification of the 9 Be bound state, and a semimicro- 
scopic continuum wave function which describes the essential scattering degrees of freedom 
at low relative energies. In Sect. 2 this model is formulated. In Sect. 3 the results for the 
9 Be(7,n) cross section are given. In Sect. 4 the astrophysical rates for a + a + n — > 9 Be+7 
and the reverse reaction are calculated. 



II. THE MODEL 

The relevant experimental data on the low-energy 9 Be(7, n) cross section are presented 
in Fig. 1. The available data are those in flTT| , |T2"f obtained with bremsstrahlung photons 
and those in |I3HI6| obtained from 7-radiation from radioactive isotopes. The peak at very 
low energy exhibited by the data of Ref. |12| is not confirmed by other groups, and may 

The radioactive 



arise from discrepancies caused by neutron energy loss in the targets [17 



isotope techniques normally provide more reliable results due to the absence of difficulties 
with the energy resolution. However, the cross section can be determined only for a few 
discrete photon energies with this method. This drawback will be cured below by the use 
of an appropriate theoretical model. Our strategy will be thus to analyze the radioactive 
isotope data. We shall consider the range of energies up to 0.5 MeV above threshold. 

We need to obtain the ground state and continuum wave functions (WF) and calculate 
the transition matrix element. We start with the three-body a + a + n representation of the 
9 Be system. Within this representation the WF in the cm. system is the aan relative motion 
function times the intrinsic WFs of the two a-particles. Since a predominant contribution 
to the transition matrix element comes from distances large compared to the size of the a- 
particle additional antisymmetrization may be disregarded.^] The intrinsic a particle WFs 
then will drop out from the calculation. In the following we shall refer to the three-particle 
relative motion function as to the WF of the system. 

Let us denote p and r the distance between the a particles and that from their center 
of mass to the neutron, respectively. The ground state is J 7 * = 3/2~, and its WF is of the 
form 



^ g s = E <PhhL(p,r)[[Y h (p)Y l2 (f)} LX i/2}j. 

hhL 



(1) 



2 The minor role of the antisymmetrization is seen from the calculations of the a + n + p system 
[18], for example. 



Here x is the neutron spin function, and the brackets [. . .] stand for angular momentum 
coupling. Because of the Pauli principle and parity requirements l\ is even, and l 2 is odd. 
The WF in Eq. (|1]) was obtained from the three-body Schrodinger equation with «a and 
an potentials reproducing the observed two-body phase shifts. These potentials, along with 
some details of the calculation, are listed in the appendix. Practically an exact solution 
to the a + a + n bound state problem is obtained, but one cannot get the experimental 
binding energy with these "bare" interparticle interactions. A possible reason is that the 
two a-particles may distort each other in the a + a + n bound state compared to the pure 
na case. This may lead to a change in the not interaction. To obtain a reasonable ground 
state WF the strength of the attractive central component of the na potential is reduced by 
8% in our calculation. This leads to values of 1.50 MeV and 2.48 fm for the binding energy 
and charge radius of 9 Be, sufficiently close to the experimental values of 1.5736 MeV and 
(2.51 ± 0.01) fm. 

Coming to the continuum wave function, we note first that the 8 Be resonance produced 
in the reaction may be safely treated as a stable particle for the purposes of our calculation. 
This is because its width of 7keV is extremely narrow on a nuclear scale. The aa continuum 
wave function, taken at the resonance energy and normalized to unity in the interior region, 
decreases practically to zero in the Coulomb barrier region as shown in Fig. 2 (solid line). 
This function represents the WF of the resonance extremely well and is taken as the 8 Be 
"bound state" wave function. Second, we argue that photodisintegration of 9 Be proceeds 
entirely into the 8 Be+n channel. Indeed, one can estimate that, at small energies consid- 
ered, the three-fragment a + a + n disintegration channel is strongly suppressed due to the 
threshold regime. The experimental data also strongly supports the absence of this channel 
TI| . At the same time, the two-fragment 4 He+ 5 He channel is still closed and ineffective 



due to the broad width of 5 He. Thus our cross section starts from the 8 Be+n threshold of 
1.6654 MeV. The cross section in the region between this threshold and the aan threshold 
of 1.5736 MeV is known to be tiny [TJ3|] and will be disregarded. As in the previous work 



(e.g. 0]) we confine ourselves to an s-wave relative n Be motion, i.e., with l/2 + continuum 
states. 

Predictions of the above dynamic aan model for the photodisintegration of 9 Be depend 
crucially on the position of the excited state of 9 Be with respect to the threshold. Preliminary 
three-body calculations gave us a peak lying too high in energy, and too broad. This could 
not be used for a reasonable fit to the data, so in the following we shall formulate an 
alternative representation of the continuum WF. We shall seek it in the form 

1 / Bc (p) 1 ^sc(r) 

*f = —ff= ~n= Xl/2- ( 2 ) 

V47T p V4vr r 

Here 8Be is the intrinsic wave function of 8 Be calculated with the same aa potential as 
for the ground state, and ip sc is the n 8 Be relative motion function, where for large r the 
normalization 

■0sc(^) — y sin(fcr + 5) (3) 

is used. 

Generally speaking, the true continuum wave function differs from Eq. (Q) not only 
in the n 8 Be interaction region but also in the outer region. However, at energies in the 
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vicinity of the long-living excited state of 9 Be the representation (Q) should approximately 
be valid in the outer region. Indeed, the decay of the long-living state into the three-body 
aan channel is inhibited due to the threshold regime. Due to the approximate validity of 
the WF in Eq. (H) in the outer region one obtains a correct energy dependence of the cross 
section when using this WF. For small energies the main energy dependence of the transition 
matrix element appears as a factor in the continuum WF and is determined by an outer 
part of the WF, i.e. the phase shift. Besides, one can see below that just the outer region 
(where tp sc takes the form of Eq. (|3|)) gives the biggest contribution to the transition matrix 
element. One can therefore hope that a WF of the form of Eq. (^) suffices for the fitting 
purpose in the whole energy range considered. 

We seek the relative motion function ip sc as a solution to the relative motion Schrodinger 
equation with some potential whose parameters are chosen from a fit of the theoretical cross 
section to the data. Taking into account that the s-wave an repulsion and the p-wave 
an attraction have comparable ranges one can assume a smooth attractive potential. The 
Woods-Saxon family 

V(r) = -Vo(l + e L ^y 1 (4) 

will be adopted below as a good representative. 

Consider now the representation of the photodisintegration cross section in our model. 
The matrix element of the dipole transition operator (4/9)er between the WFs in Eq. (|]) 
and Eq. (0) has to be calculated. After integrating over p the matrix element reduces to the 
overlap between the n 8 Be relative motion functions, namely the scattering function entering 
Eq. (H) and the "effective bound state WF" 

r- l MrWi{r)xi/2]j=z/2 (5) 

obtained from Eq. ([!]): 

r- 1 Mr)= / 0oii (p,r)prf P Be (p). (6) 



Using Eqs. (||D, (Q) the cross section is of the same form as in the single-particle case. 
The total (s-wave) photodisintegration cross section is calculated in the simple form 

a = (2/3) 6 7r(e 2 /hc)(2fi/h 2 )E^k- 1 I 2 } (7) 

/•oo 

where / = / -iph(r)^J sc (r)rdr. (8) 

•J 

Here \i is the n 8 Be reduced mass, and (hk) 2 /(2p) is the excitation energy E 1 — E th that 
will be denoted as E below. In case of a single-particle description of the process, i.e., the 
"valence neutron" model, Eq. (0) is valid with the bound-state WF normalized to unity, 
while in our case (see below) 



oo 

2 



i/j;(r)dr = 0.43. (9) 

It is implied here and in Eq. (0) that the 8 Be and 9 Be ground state wave functions are 
normalized to unity. The function ipb{ r ) is shown in Fig. 2 (dashed line). 
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III. THE PHOTODISINTEGRATION CROSS SECTION 



In Ref. [[H]] the data at the energies up to 185 keV above threshold (E y < 1.85 MeV) 



from Refs. ||13| , |14|1 were reproduced at a qualitative level within the following framework. The 
valence neutron model with Woods-Saxon type potentials was used. The results obtained 
in this way were multiplied by some constant factors less than unity, so called "reduction 
factors" , to approach the experimental cross section. To obtain the continuum wave function 
the depth of the potential was varied while the radius and diffuseness parameters were taken 
the same as for the central component of the potential in the bound state calculation. Two 
fits were found, one with a reduction factor of 0.53 leading to a weakly bound 8 Be+n l/2 + 
state, and the other one with a reduction factor of 0.31 leading to a virtual state. The 
first possibility was once preferred in view of the results of Ref. [ |22|| . In that work a two 
channel 8 Be+n model of the ground state of 9 Be was introduced to cure the single-particle 
description |20| of photodisintegration. The 8 Be subsystem was allowed to be in the ground 



and first excited state, and that led to the reduction factor of 0.5 or 0.6 in the cross section 
depending on the assumptions. In Ref. |2t| the reduction factor of 0.56 was found for that 
model. In contrast to Ref. [^1] our results below definitely testify to a virtual 8 Be+n l/2 + 



state. This is probably due to a more realistic treatment of the ground state of 9 Be in 
our model. In fact the reduction factors obtained in the two-channel model of the ground 
state of 9 Be [22| should be used in conjunction with the channel coupling n+ 8 Be dynamics, 



instead of using [^T[] single-channel dynamics. 

In Ref. p4| the same data were fitted with the line shape \/E(E + E)~ 2 . This shape was 
derived incorrectly from the Breit-Wigner cross section under the assumption that the l/2 + 



level of Be is a bound or virtual Be+n state. In Ref. [25] the data of Ref. [16| were fitted 



with a one-level R-matrix approximation. The fit leads to a complex-energy resonant state 



2q| and the real part of the complex energy proves to be negative. 

First we shall analyze the data of Refs. [p~3| — P~5f1 (full circles, full squares, and full diamond 



in Fig. 1). A search of the parameters Vo, r, and a of the potential (|4]) giving an acceptable 
fit to the data is performed. Several local minima of the quantity x 2 m the space of the 
parameters are found. One of them is provided by 

Vo = 35.99 MeV, R = 3.126 fm, a = 0.8108 fm. (10) 

These values seem to be very reasonable. For this set the X = x 2 / (degrees of freedom) value 
equals 0.62. Another one is obtained with the parameters 

V = 52.86 MeV, R = 2.006 fm, a = 1.051 fin (11) 

giving X = 0.525. Several other minima also exist with sizable higher but still acceptable X 
values. In Fig. 3 the photodisintegration cross sections obtained with the parameters (|10"D 
and (|TTJ) are shown as the solid and dashed curve, respectively. The two cross sections prove 
to be quite close to each other. To clarify partially the reason for this we note that the 
biggest contribution to the matrix element of Eq. (H) comes from the distances beyond the 
range of the potential. The distances larger than 5 fm in Eq. (||) provide 60-70% contribution 
to the cross section, and at such distances the wave functions deviate from the asymptotic 
ones, Eq. (^), by not more than 10% (except for regions in the vicinity of zeros). The 
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asymptotic wave functions are determined by the phase shifts i.e., predominantly by the 
scattering length a and the effective range r Q . Therefore the procedure is equivalent to some 
degree to fitting a and r values. Once this is done, the cross sections are not very dependent 
on the particular set of the potential parameters. The a and r values are — 27.6 fm and 
8.79 fm, respectively, for the set (|T0D, and — 28.4fm and 9.77 fm for the set (|TT]) . All the other 
above mentioned sets of potential parameters giving local minima to X lead to very similar 
a and r$ values. However considerable changes in the scattering WF inside the potential 
can influence the results, see the next paragraph. 

The following way to interpolate between the data has also been tried. Let us denote 
by cr (E) the cross section obtained in case when the right-hand side of Eq. @ is used as 
a continuum wave function for all r values. This cross section has been calculated taking 
S from the effective range expansion with the a and tq values given by the potential (pT)|). 
Let us represent a as c(E)a Q (E) and fit c(E) to experiment. It is hoped that, in contrast 
to cr, the factor c(E) behaves in a smooth way and thus can be reliably obtained from 
an interpolation procedure. Indeed, the behavior of both a and do can be approximately 
described by the resonant factor k~ l sin 2 5 times a slowly varying function. Even a fit with 
c(E) = const = 0.55 proves to provide a sufficiently low X value. The cross section obtained 
with this c is shown in Fig. 3 as the dotted curve. Presumably this procedure provides less 
accurate results than the previous one. Of course, the energy dependence k~ l sin 2 5 is not 
accurate enough in the whole energy range, as a comparison with the exact solution for the 
potential fllCf ) shows. Hence c(E) should include an energy dependence, but this could not 
be determined because of experimental uncertainties. 

Next we applied the procedure to the data of Ref. [16j (open circles in Fig. 1). Three 
local minima with acceptable X values are found. However the parameters of the potential 
corresponding to all of them: 

K) = 11.0MeV, J R = 2.35fm a = 0.258fm 

Ki = 6.49MeV, ^ = 3.10^ a = 0.260fm (12) 
V = 15.8MeV, i2=1.37fm a = 0.958fm 

prove to be rather unrealistic. There exists one more difference between these potentials 
and those in Eqs. ( [10]) and flTT|). The latter potentials, as well as the other potentials 
(|) providing a good fit to the same data, support one deeply bound state and one state 
close to being bound. On the contrary, the potentials ( JI2| ) support only one very weakly 
bound state. An existence of one deeply bound s level in the neutron mean field in the 
9 Be nucleus, or, equivalently, one node in the low-energy scattering wave function inside 
the potential, seems to be natural from the shell-model point of view. We think this point 
of view is sufficient to establish the correct number of nodes for the neutron motion inside 
the Woods-Saxon potential , even for such a clusterized system. In the a-particle oscillator 
model of 9 Be [26fl, for example, the first allowed neutron s-state contains a substantial 



admixture of the nodeless 0s function, but this leads not to a disappearance but only to a 
shift of the node. If one admits that the state considered is a mixture of 0s and Is oscillator 
functions then there exists just one node located within the distance of 3 fm from the origin. 
Therefore we conclude that in the region where various data sets differ from each other the 
older radioactive isotope data are preferable. We also note that the cross section we obtain 
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with potentials fllPf) and (|TTJ) for the highest energies considered, being lower than the fitted 
datum of Ref. ||15|| , agrees well with the bremsstrahlung Jacobson data |11 . 



IV. THE ASTROPHYSICAL REACTION RATES 

The 9 Be+7 — > a + a + n reaction rate per nucleus per time unit is calculated via the usual 
averaging the elementary photodisintegration cross section (t{E^) ■ c with the approximate, 
or Wien distribution for the photon density, 

POO 

A 7 = C 7c- 2 {hc)- 3 / a{E 1 )E 2 &qi{-E 1 /kT)dE 1 , (13) 

JE th {SBe) 

where E t h( s ~Be) = 1.6654 MeV. The rate of the reverse reaction (the number of reactions 
per time unit per unit volume) is 

P{aan) = (l/2)n 2 a n n {aan) , (14) 

where n a and n n are numbers of particles per unit volume. The reaction constant (aan) is 
obtained from Eq. (|i~3"D using the reverse ratio RR |9j27|: 



Nl(aan) = X^/RR, (15) 
RR = 5.84- 10 19 T 9 3 exp(-18.261/T 9 ). (16) 

Here N\ is Avogadro's number, T 9 is the temperature in 10 9 K, 18.261/T 9 = E t ^(aan)/kT 
with E th (aan) = 1.5736 MeV, and it is implied that the quantities ([13]) and (0) are given 
in sec -1 and cm 6 sec -1 mole -2 , respectively. Use of the Wien distribution instead of the 
exact, or Planck, one, i.e. [exp(i? 7 /fcT) — 1] _1 — > exp(— E^/kT), allows application of the 
above listed simple reverse ratio theory. For temperatures of Tg = 5 and 10, for example, 
it gives the reaction constant ( |I~3"D with relative errors of 1 % and 5.4%, respectively. For 
E 1 < 2.2 MeV the cross section o"(-E' 7 ) obtained in the preceding section with the potential 
(|I~0|) is used. For E 1 from 2.2 MeV up to 5 MeV the Jacobson bremsstrahlung data [[□] 



arc 



used. The former energy region provides 96% and 62% contribution to the cross section for 
Tg = 2 and 5, respectively. The contribution from energies E 1 higher than 5 MeV reaches 
0.3% and 6.7% for Tg = 5 and 8, respectively. 

The values of the rate (|15|) obtained can be represented by the fit 



N\(aan) = 6.59 x 10~ 6 T 9 ~ 3/2 exp[-(1.0653/T 9 )] (l + Xja„T 9 n J (17) 



-i 



with 



ai = 3.3562 a 2 = -0.86389 a 3 = 0.42268 a 4 = -0.14913 
a 5 = 2.7039 x 10~ 2 a 6 = -2.4000 x 10~ 3 a 7 = 8.3223 x 10~ 5 

The fit reproduces our N\(aan) values with the accuracy better than 1% at any Tg in the 
range 0.1 < T 9 < 8. In Eq. (|TJ) 1.0653/T 9 = [E th ( 8 Be) - E th (aan))/kT. The factor 
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Tg exp[— (1.0653/Tg)] represents the asymptotic behavior of the 8 Be formation contribu- 
tion to the rate when T 9 tends to zero. 

In the table our values for the three-body reaction rate ( |I~5"D are compared with those of 
Ref. H and those of Ref. || where constructive or destructive interference between the res- 
onant and non-resonant contributions at energies above the resonance energy was assumed. 

Summarizing, we have constructed a semi-microscopic model for the low-energy pho- 
todisintegration of the 9 Be nucleus and have analyzed the experimental data with its help. 
Our analysis supports the older radioactive isotope data. The theoretical cross section we 
derived may be compared with future microscopic calculations of the process. We have cal- 
culated the astrophysical rates for the reaction a + a + n — > 9 Be+7 and the reverse reaction. 
Our new reaction rates agree at T 9 = 2.0 with the ones given in Ref. |J. They are somewhat 
smaller (larger) for lower (higher) temperatures than T 9 = 2.0. The reaction rates given 
in Ref. || agree much better with our reaction rate at higher temperatures if one assumes 
in Ref. |J constructive (destructive) interference between the resonant and non-resonant 
contributions at energies above (below) the resonance energy. 

We are indebted to J.S. Vaagen and J.M. Bang for very fruitful comments. This work 
was supported partially by the the Fonds zur Forderung wissenschaftlichen Forschung in 
Osterreich (project P10361-PHY) and the Russian Foundation for Basic Research (grant no 
97-02-17003). 



APPENDIX 

In our calculation of the ground state of 9 Be the aa potential is taken in the form 

V R exp[-(p R p) 2 ] - V A ex\o[-(/j, A p) 2 ] 

with V A = 130 MeV, p A = 0.475 fm" 1 , p R = 0.7 fin" 1 , and V R = 500, 320, and lOMeV for 
1=0, 2, and 4, respectively. The Coulomb aa interaction is also added. The not- interaction 
in s-, p-, and d-states is taken into account. As in many previous studies BO] the s-wave 



repulsive potential Vexp[— (r/R) 2 ] with V = 50 MeV and R = 2.3 fm is used. The initial 
potential in p- and d-states |29| includes central and spin-orbit components: 



/ r-R\~l Id/ r - R l 

V(r) = —V (1 + e~ ) - U 1 • s -— 1 + e »i 
v ' r dr V 

with V = 43 MeV, R = 2fm, a = 0.7 fm, U = 40MeVfm 2 , R x = 1.5 fm, and a x =0.35 fm. 
The parameter V is reduced to 39.6 MeV in the present three-body calculation, in order to 
reproduce the empirical g.s. energy. 

The three-body dynamic equation is written in the form of the Faddeev differential 
equations and each Faddeev component is expanded over hyperspherical harmonics and 
hyperradial basis functions. Using the Raynal-Revai rotations of hyperspherical harmonics^ 
the matrix elements are reduced analytically to two-dimensional integrals. The equations 



3 A review of the hyperspherical formalism can be found e.g. in Ref. puj\ . 

9 



are projected onto subspaces of the basis functions retained that reduces the problem to 
the algebraic eigenvalue problem. The number of basis functions retained is quite high and 
ensures the adequate convergence of the calculation. 



10 



REFERENCES 



[1] A.S. Pudliner, V.R. Pandharipande, J. Carlson, S.C. Pieper, and R.B. Wiringa, Phys. 

Rev. C56, 1720 (1997). 
[2] K. Arai, Y. Ogawa, Y. Suzuki, and K. Varga, Phys. Rev. C54, 132 (1996). 
[3] S.E. Woosely, J.R. Wilson, G.J. Mathews, R.D. Hoffmann, and B.S. Meyer, Astrophys. 

J. 433, 229 (1994). 

[4] K. Takahashi, J. Witti, and H.T. Janka, Astron. and Astrophys. 286, 857 (1994). 
[5] S.E. Woosely and R.D. Hoffman, Astrophys. J. 395, 202 (1992). 

[6] J. G6rres, H. Herndl, I.J. Thompson, and M. Wiescher, Phys. Rev. C 52, 2231 (1995). 
[7] V.D. Efros, W. Balogh, H. Herndl, R. Hofinger, and H. Oberhummer, Z. Phys. A 355, 
101 (1996). 

[8] H. Herndl, R. Hofinger, and H. Oberhummer, in Editors: S. Kubono and T. Ka- 
jino, Proceedings of the Symposium on Origin of Matter and Evolution of Galaxies 
97 (OMEG97), Nov. 5-7, 1997, Atami, Japan, World Scientific, Singapore, in press. 

[9] W.A. Fowler, G.R. Caughlan, and B.A. Zimmerman, Ann. Rev. Astron. Astrophys. 13, 
69 (1975). 

[10] F. Ajzenberg-Selove, Nucl. Phys. A490, 1 (1988). 
[11] M.J. Jacobson, Phys. Rev. 123, 229 (1961). 

[12] B.L. Berman, R.L. Van Hemert, and CD. Bowman, Phys. Rev. 163, 958 (1967). 
[13] J.H. Gibbons, R.L. Maclin, J.B. Marion, and HW. Schmitt, Phys. Rev. 114, 1319 
(1959). 

[14] W. John and J.M. Prosser, Phys. Rev. 127, 231 (1962). 

[15] B. Hammermesh and C. Kimball, Phys. Rev. 90, 1063 (1953). 

[16] M. Fuishiro, T. Tabata, K. Okamoto, and T. Tsuimoto, Can. J. Phys. 60, 1672 (1982). 
[17] F.C. Barker and B.M. Fitzpatrick, Aust. J. Phys. 21 415 (1968). 

[18] G.G. Ryzhih, R.A. Eramzhyan, V.I. Kukulin, and Yu.M. Tchuvilsky, Nucl. Phys. A563, 
247 (1993). 

[19] M. Fuishiro, K. Okamoto, and T. Tsuimoto, Can. J. Phys. 61, 1579 (1983). 
[20] N.C. Francis, D.T. Goldman, and E. Guth, Phys. Rev. 120, 2175 (1960). 
[21] E.G. Corman, J.E. Sherwood, and W. John, Phys. Lett. 4, 146 (1963). 
[22] J.S. Blair, Phys. Rev. 123, 2151 (1961). 
[23] F.C. Barker, Nucl. Phys. 28, 96 (1961). 
[24] C. Mahaux, Nucl. Phys. 71, 241 (1965). 
[25] F.C. Barker, Can. J. Phys. 61, 1371 (1983). 
[26] PD. Kunz, Ann. Phys. 11, 275 (1960). 

[27] W.A. Fowler, G.R. Caughlan, and B.A. Zimmerman, Ann. Rev. Astron. Astrophys. 13, 
69 (1975). 

[28] S. Ali and A.R. Bodmer, Nucl. Phys. 80, 99 (1966). 
[29] J. Bang and C. Gignoux, Nucl. Phys. A313, 119 (1979). 

[30] M.V. Zhukov, B.V. Danilin, D.V. Fedorov, J.M. Bang, I.J. Thompson, and J.S. Vaagen, 
Phys. Rep. 231, 151 (1993). 



11 



TABLES 



1 [10 KJ 






JN\{aan) 






Present work 


Ref. |9J 


Ker. |gj (destructive) 


Ket. |qj (constructive) 


0.2 


0.21-10- 6 


0.30-10" 6 






0.5 


0.87-10- 6 


1.1-10 -6 


0.55 • 10~ 6 


0.61 • 10~ 6 


1.0 


0.60-10- 6 


0.67-10- 6 


0.32 • 10~ 6 


0.44 • 10~ 6 


2.0 


0.23-10" 6 


0.23-10~ 6 


0.12 • 10~ 6 


0.20 • 10~ 6 


3.0 


0.12-10" 6 


0.99-10- 7 


0.60 • 10~ 7 


0.11 • 10~ 6 


4.0 


0.73-10- 7 


0.52-10- 7 






5.0 


0.51-10- 7 


0.31-10- 7 


0.26 • 10~ 7 


0.52 • 10~ 7 



12 



FIGURES 



FIG. 1. Experimental data on the low-energy cross section for photodisintegration of 9 Be. 
Bremsstrahlung data: Ref. [12] (solid curve representing the fit of the authors to their data), 
Ref. [11] (star representing the maximum of the spectrum; the rest of the spectrum is not shown). 
Radioactive isotope data: Ref. [13] (full circles), Ref. [14] (full squares), Ref. [15] (full diamond), 
Ref. [16] (open circles) 

FIG. 2. The wave function of the 8 Be resonance (solid curve) and the "effective" n- 8 Be relative 
motion wave function in 9 Be (dashed curve). The former wave function is constructed as the 
s-wave continuum solution in the aa potential (see appendix) at the energy of 0.09518 MeV, the 
peak of the resonance in this potential. 

FIG. 3. Calculated cross sections for the photodisintegration of 9 Be. Solid and dashed curves 
are obtained with the parameters (10) and (11), respectively. Dotted curve corresponds to the 
alternative method listed in the paragraph after Eq. (11). The experimental data are as in Fig. 1. 
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